Shedding light on the calculation of electrode electroactive area and heterogeneous electron transfer rate constants at graphite screen-printed electrodes

We present in detail the most known and commonly used methods for the calculation of electrode electroactive area (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document}A) and heterogeneous electron transfer rate constants (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}^{0}$$\end{document}k0). The correct procedure for the calculation of these parameters is often disregarded due to either lack of a minimum theoretical background or oversimplification of each method’s limitations and prerequisites. The aim of this work is to provide the theoretical background as well as a detailed guide for the implementation of these measurements by impressing upon the electrochemists the parameters that need to be considered so that the obtained results are safe and useful. Using graphite screen-printed electrodes, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document}A, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}^{0}$$\end{document}k0 were calculated with different methods and techniques. Data are compared and discussed. Graphical Abstract


Calculation of electrode electroactive area (A)
An important parameter to consider when evaluating the response of an electrochemical system is the electroactive area of the working electrode as the peak current ( I p ) is dependent on A regardless of the system's behavior (reversible, quasi-reversible or irreversible). The A of an electrode is also correlated with the current density j (A cm 2 ) observed during an electrochemical experiment, and thus its knowledge is necessary for obtaining the full picture of the electrochemical process under investigation. It is also involved in most calculations of other electrochemical parameters ( k 0 , C dl , etc.) and properties (for example, the electrocatalytic performance of a material) since most of them are often (and should be) normalized with A [1][2][3][4]. Setting a clear pathway for its calculation is of paramount importance as it is probably one of the first studies an electrochemist should embark upon when evaluating a new system. Finally, as regards the screen-printed electrodes (SPEs) on which the case study is focused, it should be noted, that even though the interelectrode reproducibility of a single batch of SPEs is rather satisfactory, it is highly recommended that the calculation of A is performed at every new batch, as the status of both the ink and of the printing mesh as well as the printing settings can influence its value. This also applies for other types of non-conventional electrodes, such as the 3D-printed and the laser scribed electrodes, the A of which is greatly affected by the fabrication settings, the source material etc.
There are two main techniques for the calculation of the electroactive area: i) chronocoulometry and ii) cyclic voltammetry, which will be presented in detail.

Chronocoulometry
Chronocoulometry is a simple and commonly employed technique for the calculation of the electroactive area of Maria G. Trachioti and Alexandros Ch. Lazanas contributed equally to this work. an electrode, by performing experiments in the presence of a redox compound of known diffusion coefficient ( D ). As its name implies, chronocoulometry is the measurement of charge ( Q , in coulombs) with time (chrono) and is implemented by applying a short dc potential step to an electrochemical cell, while the current (I, in amperes) transient is monitored. Charge versus time plots are provided by the electrochemical analyzers by integrating the current over time. In double potential step chronocoulometry, the potential of the working electrode (in fact, the potential difference versus the reference electrode) is stepped from an initial value ( E i ) where no electrochemical Faradaic reaction occurs to a value ( E s ) where complete electrolysis of the redox species occurs, so that the current is limited by (planar) diffusion, and then to a final value ( E f ), that is frequently identical to E i [5][6][7]. Obviously, experiments are conducted under quiescent conditions. The transition from E i to E s is commonly referred to as forward step, and from E s to E f (= E i ), as reverse step. The duration of each potential step is taken as , in seconds.
The Anson equation (Eq. 1) calculates the total charge (Q T ) , which passes through the electrochemical cell at any time, due to the charging of the electrical double layer (Q dl ) and the Faradaic reaction of diffusing ( Q diff ) and/or absorbed ( Q ads ) species of the redox compound, and can be applied to the calculation of the electroactive area of the working electrode (A) [6,8,9].
Redox species are classified as outer-or inner-sphere depending on the way the electron transfer processes occur on the electrode's surface [10]. In the case of outer-sphere redox species, the electrochemical process is diffusion controlled and thus is influenced by the electronic properties of the electrode surface only, while in the case of inner-sphere redox species, the electrochemical process is also influenced by surface functional groups and commonly involves adsorbed species [10][11][12].
In this regard, the Anson equation for a potential step when both diffusing and absorbed species are included in the electrochemical process is given by the following equation where Q T is the total charge (in coulombs, C ), n is the number of electrons involved in the electrochemical reaction, F is Faraday's constant (96485 C moL −1 ), A is the electrode's electroactive area in cm 2 , C is the concentration of the redox species in solution (moL cm −3 ), D is the diffusion coefficient (cm 2 s −1 ), while Q dl , Q ads and Q diff are the charge components due to charging of the electrical double layer, the electrolysis of adsorbed species, and the electrolysis of solution (diffusing) species, respectively. The charge component Q ads , if any, is equal to nFAΓ 0 , where Γ 0 is the amount of adsorbed species on the electrode surface.
To discriminate among the three components of the total charge, Q dl , Q ads and Q diff , a plot of Q T versus √ t (during the forward step) and ( √ + √ t − − √ t) (during the reverse step), collectively termed as Anson plot, can be employed. Considering that Q dl and Q ads remain constant with the time, by using the Anson plot's slope attributed to the diffusioncontrolled charge transfer process (for an ideally reversible process the slopes at both forward and reverse steps should be the same), the electroactive area of the electrode can be calculated as per the following equation

Cyclic voltammetry
Considering the factors that should be taken into account regarding the size and the morphology of the electrodes (for a detailed analysis of these factors the reader is referred to the Refs [11,[13][14][15]) if cyclic voltammetry is employed, the electrode electroactive area can be calculated, in the presence of a redox species in the measuring solution, by using the proper variation of the Randles-Ševčík equation.
In the case of a reversible process ( nΔE p ≈ 57.5 mV, ΔE p is the peak-to-peak potential separation), the peak current is described by the Randles-Ševčík equation where I p is the forward peak current, v is the potential scan rate, while the other terms have the aforesaid meaning. Thus, for a reversible system, the peak current is proportional to the square root of the scan rate. A plot of I p versus √ v gives a straight line, the slope of which can be used to calculate the A [2,5,9].
Nevertheless, despite the popular bibliographical belief, common redox molecules (for example, potassium ferrocyanide) are more often involved in quasi-reversible processes (typically, 63 < nΔE p < 200 mV), especially when non-conventional electrode materials, such as the widely used in modern electro-sensing applications, screen-printed or 3D-printed electrodes, are employed. In these cases, if 63 < nΔE p < 200 mV, Eq. 3 is not valid and the modified Randles-Ševčík equation for quasi-reversible processes should be used instead [16] (2) where K(Λ, ) is a modified dimensionless parameter for quasi-reversible reactions. To estimate the K( , ) parameter, we have first to calculate another dimensionless parameter termed ( , ) by using Eq. 5, and parameter by using Eq. 6 [17] where E p is the peak potential, E p∕2 is the half-wave potential (these values can be found from the recorded cyclic voltammograms), R is the gas constant (8.314 J mol −1 K −1 ) and T is the temperature ( K).
where Ψ is a kinetic parameter and its calculation is described analytically below. The next step is a simple calculation of logΛ . Using the plot Δ(Λ, ) = f (logΛ) (Fig. 1A), as the parameters Δ(Λ, ) and logΛ are known, the transfer coefficient ( ) can be approximated graphically. Finally, the parameter K(Λ, ) is found graphically from the plot K(Λ, ) = f (logΛ) for the estimated value of ( Fig. 1B) [2,18]. Note: We will see below (eqs. 11 and 12) that the Ψ parameter is given as a function of ΔE p , which in turn, for a non-reversible process is dependent on the scan rate. Since value emerges for each different scan rate value [2]. In other words, the K(Λ, ) parameter in Eq. 4, declares that at a quasi-reversible process ( 63 < nΔE p < 200 mV), the peak current is not proportional to the square root of scan rate.
In the case of an irreversible process (typically, when nΔE p > 200 mV), the peak current is given by the following equation where n and n ′ are the number of electrons in the electrochemical reaction and the number of electrons transferred before the rate determining step, respectively [11,19], while all the other terms have their aforementioned meaning. Assuming that n = n � = 1 , the peak current is calculated by using the modified Randles-Ševčík equation given below The transfer coefficient ( ) is a measure of symmetry of the energy barrier for a single electron transfer step [2]. For a purely symmetrical reaction in terms of energy (a reversible reaction) has a transfer coefficient of 0.5 [1,19]. Consequently, in irreversible reactions should not be considered to be 0.5. In this case, can be calculated by the Eq. 9 [2]

Calculation of heterogeneous electron transfer rate constants (k 0 )
Heterogeneous electron transfer rate constants ( k 0 ) have been a key parameter of the electrochemical performance of various electrode materials, or electrode modifications, or electrode modifiers since they reflect the kinetics of the reaction between a particular redox compound and the electrode surface [19]. In general, k 0 is a measure of the heterogeneous kinetic facility between the oxidized and the reduced forms of a redox couple. An electrochemical system with a large k 0 can achieve equilibrium quickly, while a system with small k 0 , sluggishly [2]. Admittedly, cyclic voltammetry is by far  [21], Gileadi [22], etc.) have been developed on it. Besides, even though not so widely used, electrochemical impedance spectroscopy (EIS) can also be used for the calculation of k 0 both in a numerical [16,[23][24][25] and in a graphical fashion [26]. Note: The rate constant of an electrode reaction refers to the time needed for the electroactive species to arrange themselves and their ionic atmospheres for the electron transfer to occur. It does not measure the rate of electron transfer itself, as this occurs extremely rapidly, in approximately 10 −16 s [1].

Nicholson method [20]
One of the most commonly used voltammetry-based methods is Nicholson's method [20], in which k 0 for quasi-reversible electrochemical reactions is related to a dimensionless kinetic parameter, named Ψ , as per Eq. 10 where Ψ works as a polynomial function of the product nΔE p within the range 63 < nΔE p < 212 mV, or with respect to Ψ values, within the ( 7 > Ψ > 0.1 ) range. The upper Ψ limit of 7 (or nΔE p > 63 mV) signifies the transition between a reversible to a quasi-reversible system, while the lower limit of 0.1 (or nΔE p > 212 mV) signifies the transition between a quasi-reversible system to an irreversible one. Practically, the only input required to calculate k 0 is the value of the Ψ parameter, provided that the nΔE p value is within the quasireversible limits ( 63 < nΔE p < 212 mV) set above. Obviously, the nΔE p value is highly dependent on the experimental variables (the concentration of the redox compound, the electrolyte, the pH of the measuring solution, the scan rate, etc.), which can thus be appropriately tuned to set nΔE p within the quasi-reversible limits ( 63 < nΔE p < 212 mV), where the Nicholson method is valid.
The Ψ parameter can be found graphically from Fig. 2 (this graph represents a better resolution graph of Fig. 3 in [20] and it was constructed based on the data given in Table 1 of [20]). More conveniently, Ψ parameter can be calculated by using the empirical equation developed by Lavagnini et al. [27]. In fact, Eq. 11 is a polynomial fit of Nicholson's working curve and allows the easy and accurate determination of Ψ parameter from the nΔE p of a single cyclic voltammogram. When nΔE p ranges between 140 and 200 mV, the same authors [27] suggest that Ψ parameter to be calculated by using the following equation which offers a more representative Ψ value in this range. Note: An alternative graphical approach is plotting the Eq. 10 as an f (x) = bx function, where f (x) = Ψ , x = 1∕ √ Dn F∕RT and b = k 0 (slope) at different scan rate values, with the restriction of nΔE p being kept under 212 mV. As a general rule, the average value of k o at different scan rates is considered more reliable than its value at a single scan rate.

Klingler-Kochi method [21]
In the case of irreversible systems ( nΔE p > 150 mV), the following equation suggested by Klingler and Kochi [21] is available and should be considered: where , the transfer coefficient of the forward scan, is calculated as described above (Eq. 9) and all the other symbols have the same meaning.
Note: when nΔE p > 150 mV, by using the Klingler and Kochi method [21] and Eq. 13, k 0 can be directly calculated from a single experimental variable ( nΔE p )  Table 1 in Ref. [20] Page 5 of 14 251 without being necessary any other parameter (for example, Ψ ) to be previously calculated. Remember though that the average value at different scan rates is preferred.

Gileadi method [22]
In this method, there is no restriction to the ΔE p value since it is based on the graphical approximation of the critical scan rate (at which the electrochemical process changes from reversible or quasi-reversible to irreversible) from a plot presenting the variation of the peak potential with the logarithm of scan rate [ E p = f (logv) ]. Two straight lines with different slopes are obtained at the low and high scan rates and the value of the critical scan rate is estimated from the intersection of the extrapolated lines. Then, the calculation of k 0 is performed using the following equation where v c is the critical scan rate and the other terms have their aforesaid meaning [22,28]. Note: The methods presented above, and especially the Gileadi method which is based on the graphical calculation of k 0 should not be confused with the well-known method of Laviron [29], which refers to surface confined redox species and as a result the respective rate constant is given in s −1 and not in cm s −1 .

EIS-based methods
EIS enables the calculation of k 0 through the determination of charge transfer resistance ( R ct ) when EIS measurements are conducted in the presence of a redox couple (the electrolyte contains both the Ox and Red forms of a redox system).
The relationship between k 0 and EIS derived data has been introduced by Randles [30], and has been studied extensively by Sluyters [23]. The apparent limitation of this method is that the impedance spectrum must be sufficiently modeled by a Randles equivalent electrical circuit [26] and that the overpotential of the redox reaction should be sufficiently low (see Results & Discussion). Note: k 0 may be also determined by EIS when examining complex electrochemical that cannot be simulated to a single Randles circuit. In those cases, profound knowledge of the electrochemical system at hand is essential because the existence of additional charge-transfer phenomena can lead to the existence of more than one R ct values in the equivalent circuit.
The estimation of k 0 by EIS data is possible when working at the linear part of the Butler-Volmer equation, valid for small overpotential values [2]. Under these conditions, the exchange current i 0 is related with the (small) values of R ct with the Eq. 15 Exchange current i 0 is also correlated with k 0 according to Eq. 16: where C * = C * ox = C * red , is the concentration of the redox couple in the solution. By combining eqs. 15 and 16, the correlation between R ct and k 0 is enabled as follows: at increasing concentrations of the redox probe, which will result at decreasing R ct [26].
Note: For an ideally reversible reaction the overpotential is zero and consequently so is the exchange current. By extension, that means that the R ct is also zero and k 0 tends to infinity ( R ct =0, k 0 → ∞ ) [8].
An overview of the various techniques and methods elaborated in this case study is summarized in Table 1, while the experiments conducted, and the respective results are presented in Table 2.

Materials
Hexaammineruthenium(III) chloride (RuHex) was purchased from Aldrich. Potassium hexacyanoferrate(III) (ferricyanide) was purchased by AnalaR. Potassium hexacyanoferrate(II) trihydrate (ferrocyanide) was a product of Merck. All the other reagents were of analytical grade from Sigma-Aldrich.

Apparatus
Electrochemical measurements were conducted with a PGSTAT12/FRAII electrochemical analyzer (Metrohm and not by applying the equation to a single scan rate is suggested Quasi-reversible process: 63 < nΔE p < 200 mV K(Λ, ) is found as follows: i) calculation of Ψ (see below in this column), ii) calculation of logΛ as Λ = Ψ √ , iii) calculation of Δ(Λ, ) as Δ(Λ, ) = E p∕2 − E p ∕26, iv) a can be found graphically from the plot Δ(Λ, ) = f (logΛ), v) K(Λ, ) can be found graphically from the plot K(Λ, ) = f (logΛ) for the a value approximated in step iv Irreversible process: v c is estimated graphically through the E p = f (logv) plot Electrochemical Impedance Spectroscopy k o = RT∕n 2 F 2 R ct AC On condition that impedance spectrum fitting Randles circuit and that the overpotential of the redox probe is low (typically < 50 mV; see

Fabrication of graphite SPEs
The graphite SPEs (4 mm diameter) were fabricated in arrays of twenty electrodes onto a 175 μm polyester substrate (Autostat CUS7) using a semi-automatic screen printer (E2, EKRA) and home-made polyester screens with 195 mesh (77/195-48 PW, SEFAR PET 1500). Graphite ink (Loctite EDAG 407A) was pushed into the open area of the screens with a 75-durometer polyurethane squeegee. SPEs were cured at 90 °C for 60 min in a conventional oven. Table 2 An overview of the techniques and methods used in this case study, as well as the obtained results. Values are presented as mean ± SD, n = 3. CC, chronocoulometry; CV, cyclic voltammetry; EIS, electrochemical impedance spectroscopy

Calculation of electroactive area (A)
Chronocoulometric and cyclic voltammetry experiments using two different redox probes (an outer-sphere, RuHex, and an inner-sphere, ferricyanide) were conducted, and the electrode electroactive area was calculated with the various techniques/methods discussed above. Calculated values are also expressed as the percentage ratio of the calculated area ( A ) to the geometric area ( A geo ) ( %R = A∕A geo × 100 ) [11]. The A geo of SPEs (d = 0.4 cm) is 0.1256 cm 2 .
As regards the electroactive area calculated using double potential step chronocoulometry, measurements were conducted in 1 mM RuHex in 0.1 M KCl or 1 mM ferricyanide in 0.1 M KCl, pH 3 or 0.1 M PBS, pH 7. The values of E s and E i = E f for RuHex were − 0.5 V and 0.2 V and for ferricyanide − 0.4 V and 0.8 V, respectively. The total charge that passed through the cell was recorded for 6 s (3 s at E s and 3 s at E f ) to ensure that the redox probe will diffuse to the electrode at its maximum rate [2], in accordance with previous studies [11]. Chronocoulometric measurements (3 s at E s and 3 s at E f ) were preceded by a conditioning step at E i for 3 s. The experiments were conducted with 3 different SPEs. The signal recorded in the absence of the redox probe (in pure electrolyte) was subtracted from the signal recorded in the presence of the redox probe. The subtracted chronocoulometric signals of three different SPEs in 0.1 M KCl containing 1 mM RuHex and the respective Anson plots are presented in Fig. 3. Generally, while the calculation of A is feasible using either slope of the Anson plot, we have selected the one that corresponds to the E f in each case, since it presents the most reproducible results, and the linearity of the respective curves is better. The diffusion coefficients used in this work were 9.1 × 10 −6 and 7.6 × 10 −6 cm 2 s −1 for RuHex and ferricyanide, respectively. The electroactive area calculated using RuHex is 0.1154 ± 0.0032 cm 2 ( %R = 91.9%), while the areas calculated using ferricyanide in KCl, pH 3 and in PBS, pH 7 are 0.1510 ± 0.0052 ( %R = 120.2%) and 0.1481 ± 0.0087 cm 2 ( %R = 117.9%), respectively.
To calculate the electrochemically active area of the electrodes via cyclic voltammetry, by using, depending on the ΔE p for each recorded CV, the appropriate in each case Randles Fig. 4, while to compare the CV responses at the three examined measuring solutions, the CVs obtained in each case at the same scan rate value (50 mV s −1 ) are presented in Fig. 5. From the CVs in 0.1 M KCl containing 1 mM RuHex, the average electroactive area of 3 SPEs was calculated for scan rates 5 − 300 mV s −1 , where 63 < ΔE p < 200 mV, using the Randles-Ševčík equation for quasi-reversible processes (Eq. 4) and was found 0.1002 ± 0.0023 cm 2 ( %R = 79.8%). From the CVs in ferricyanide in 0.1 M KCl, pH 3, average A was calculated from the scans conducted at scan rate 5 − 50 mV s −1 , with the modified Randles-Ševčík equation for quasi-reversible processes (Eq. 4), while from the CVs conducted at scan rates 75 − 500 mV s −1 , due to the large ΔE p at the respective CVs, with the Randles-Ševčík equation for irreversible processes (Eq. 8). A was calculated 0.0990 ± 0.0023 ( %R = 78.8%) and 0.0941 ± 0.0058 cm 2 ( %R = 74.9%), respectively. From the measurements in ferricyanide in 0.1 M PBS pH 7, due to the large ΔE p at all the examined scan rates, average A was also calculated by using the Randles-Ševčík equation for irreversible processes (Eq. 8) and was found to be 0.0331 ± 0.0035 cm 2 ( %R = 26.3%). a was calculated for each scan rate value at which the cyclic voltammetric response become irreversible (at high scan rates in ferricyanide in 0.1 M KCl, pH 3 or at all scan rates in ferricyanide in 0.1 M PBS, pH 7) using Eq. 9 and was used accordingly.
Concerning the values of A calculated by chronocoulometric experiments with RuHex in KCl and with ferricyanide in 0.1 M KCl, pH 3 or 0.1 M PBS, pH 7, the former was found to be near to A geo ( %R = 91.9%) in contrast with the two values obtained with ferricyanide, which are considerably higher ( %R = 120.2 and 117.9%).
From the values of A calculated by cyclic voltammetry measurements, the experiment conducted with ferricyanide in 0.1 M PBS, pH 7 seems to be a poor choice ( %R was only 26.3%), while the results obtained from the other experiments could be considered satisfying. Data indicate that the modified Randles-Ševčík equation for irreversible systems cannot provide safe results for the calculation of A , especially when the ΔE p values are too high (500 to 800 mV). This is not the case, however, for the ferricyanide in 0.1 M KCl, pH 3, especially at low scan rates (5-50 mV s −1 ) where the respective ΔE p values indicate a quasi-reversible response.
To summarize, the electroactive area calculated with RuHex via chronocoulometry, and the respective values calculated with the same redox probe via cyclic voltammetry, as well as with ferricyanide in 0.1 M KCl, pH 3 at low scan rates are quite similar. Even though we cannot compare the obtained results with respect to a reference value, chronocoulometric measurements using RuHex seems to be more advantageous than the others for the calculation of the electroactive area of graphite SPEs, since it offers fast and reliable results (considering the calculated %R in each case), in agreement with previous works [11].
Quite often in the literature, in the case of quasi-reversible processes the calculation of A has been conducted incorrectly using the Randles-Ševčík equation for reversible processes, and not by using the modified Randles-Ševčík equations. This predication that the I p is proportional to the √ v in any case leads to an important error. Additionally, in irreversible reactions a similar error occurs when is considered 0.5, and not calculated by the E p and E p∕2 values (see Eq. 9) at each examined scan rate. To clear these resulting errors up, we demonstrate the use of a I p ∕L = f � √ v � function, where " L " (which in this case L stands for linearity) is the necessary divisor for maintaining the linear dependence of I p to the √ v (Fig. 6). More specifically, in the case of quasi-reversible processes the L factor is equal to the K(Λ, ) (see Eq. 4) and in the case of irreversible processes the factor L is equal to 1.11 √ (see Eq. 8). The discrimination of the electrochemical process as "quasi-reversible" or "irreversible" is based on the value of nΔE p at each CV. When nΔE p < 200mV , the I p value is corrected to the corresponding K(Λ, ) value as dictated by Eq. 4 for a quasi-reversible process, and when nΔE p > 200mV , the I p value was corrected to 1.11 √ (Eq. 8) while the value of was calculated according to the Eq. 9 from the E p and E p∕2 values at each scan rate. Obviously, in reversible processes the factor L is equal to 1 . As it is evident from the Fig. 6, correcting each experimental value of I p with the respective L factor is indispensable for the accurate calculation of A . Thus, the calculation of A from the slope of this linear plot, over all the examined scan rates, both for the RuHex and for ferricyanide in 0.1. M KCl, pH 3, was conducted and the respective values (0.1015 ± 0.0010 and 0.0951 ± 0.0021 cm 2 ) are almost the same with the values calculated from the respective (depending on the nΔE p at a given scan rate) Randles-Ševčík equations.

Heterogeneous electron transfer rate constants (k 0 )
The same CV experiments used for the calculation of electrode electroactive area with the two different redox probes (RuHex and ferricyanide) can also be exploited for the calculation of heterogeneous electron transfer rate constants.
Regarding the value of ΔE p in each cyclic voltammogram at a specific scan rate, the value of k 0 was calculated either with the Nicholson method ( ΔE p < 200 mV) [20] or with the Klingler-Kochi method ( ΔE p > 150 mV) [21]. Ψ values (used in Nicholson method) were calculated by Eq. 11 and Eq. 12, for ΔE p < 150 mV and ΔE p > 150 mV, respectively. Specifically, for RuHex the values of k 0 were found (2.535 ± 0.039) × 10 −3 (Nicholson method), for scan rates 5 -300 mV s The difference in the k 0 values calculated by the Nicholson method and the Klingler-Kochi method for the case of ferricyanide in 0.1 M KCl, pH 3 can be explained by the fact that the Nicholson method is based on the approximation of the kinetic parameter Ψ for =0.5, without taking into consideration the true value of , especially at high scan rates. Indeed, as stated in the original work by Nicholson, as Ψ becomes sufficiently low (for large ΔE p values), the error in the calculation of Ψ due to the different values of (which in this case should not be considered to be 0.5) spans from 5% ( Ψ =0.5 or ΔE p = 105mV) to 20% ( Ψ = 0.1 or ΔE p = 212mV) . On the other hand, the Klingler-Kochi method is based on the fact that the electrochemical system should be forced to the irreversible regime where the influence of is sound and thus requires the experimental value of at a given scan rate and considers this value in Eq. 13. For this reason, for nΔE p > 150 mV the use of Klingler-Kochi method is more suitable for the calculation of k 0 (as seen in Eq. 13, the value of is considered), while the Nicholson method remains the best choice for nΔE p values up to ca. 140 mV (ideally up to 105 mV) where can be considered to play a minimum role in the kinetics of the system.
Here it is important to note that contrary to the case of ferricyanide in 0.1 M KCl, pH 3, the two methods gave almost identical results in the case of RuHex in 0.1 M KCl. This is due to the fact that the true values of , as they were found experimentally at each scan rate, were very close to the theoretical value of 0.5.
Then, the value of k 0 was estimated with the graphical Gileadi method [22]. This method is not applicable in the case of ferricyanide in 0.1 M PBS, pH 7, because the process is irreversible at every applied scan rate. On the other hand, with the other two examined systems, the estimation of k 0 with the Gileadi method is feasible and the values of k 0 were found to be (3.509 ± 0.143) × 10 −3 (RuHex) and (3.310 ± 0.250) × 10 −3 cm s −1 (ferricyanide in KCl, pH 3), respectively. At these two cases, due to the transition of the electrochemical process from quasi-reversible to irreversible as the scan rate increases, it is easy to estimate the value of v c (i.e., the value where this transition happens). With the same reasoning as for the Klingler-Kochi method, we can deduct that the Gileadi method produces values of k 0 close to that of Klingler-Kochi method because it considers the variation of (see Eq. 14) both when approaching and after progressing to the irreversible regime. An indicative graph constructed for Gileadi method in the case of RuHex is depicted in Fig. 7.
The estimation of k 0 using EIS data was conducted, as well. The aforementioned ferro/ferricyanide couple in 0.1 M KCl, pH 3 was not the optimum choice due to the system's instability in successive EIS measurements, as can be seen in Fig. 8A. According to previous studies, this behavior can be attributed to specific adsorption phenomena [31][32][33][34]. On the other hand, in 0.1 M PBS pH 7, EIS spectra are quite stable over time. Thus, EIS measurements were conducted in a mixture of 1 + 1 mM ferro/ferricyanide in 0.1 M PBS, pH 7 and the respective spectra (three successive scans) are shown in Fig. 8B. The EIS excitation signal was superimposed to a dc potential of 0.210 V (equal to formal potential E 0 = E pa + E pc ∕2 , which was found by conducting a CV at 50 mV s −1 ). The composition of the measuring solution and the DC potential were selected to ensure equimolar concentrations of the oxidized and reduced species at the electrode surface at all times ( C * = C * ox = C * red ) . EIS data were fitted to a Randles circuit (Fig. 8B, inset), which was used to estimate the value of R ct (Table 3) from which (using the Eq. 17) k 0 was estimated (0.022 ± 0.002) × 10 −3 cm s −1 .
This value is an order of magnitude lower than the respective value found by the CV-based method for the same system ( Table 2). The explanation for this deviation lies on the fact that Eq. 17 is valid for very low overpotential ( ) values as a result of a linear approximation of the Butler-Volmer equation. If no mass-transfer effects are considered, the Butler-Volmer equation is [2]: while the approximated linear version of the Butler-Volmer equation is [2]: By definition, is the additional voltage required in a non-reversible electrochemical system (by reference to its formal potential when C * ox = C * red ) for the electrochemical reaction to occur. For example, for the anodic part of a redox reaction is: It is evident that in this specific system (ferro/ferricyanide in 0.1 M PBS, pH 7) the overpotential value is high enough (> 200 mV) to negate the use of Eq. 17 for the calculation of k 0 via EIS measurements since in this case, R ct and k 0 are not linearly related. This begs the question about the maximum value that permits the use of Eq. 17 and consequently, the calculation of k 0 via EIS data in a particular electrochemical system (redox molecule, electrolyte, electrode, etc.). In response, we provide a graph illustrating the current values generated   (Fig. 9). The highlighted areas show the average error in current, I error (%) = [(I B.V. − I L )∕I B.V. ] × 100 , and consequently the validity of the linear approximation with respect to . Within the green area ( spans from 0 to ± 50 mV) the mean error is 5%, within the yellow area ( spans from ± 50 to ± 100 mV) 31%, within the orange area ( spans from ± 100 to ± 200 mV) 69%, and finally, within the salmon area ( > ±200 mV) > 86%. With that in mind, we can deduct that a suitable electrochemical system to use for the estimation of k 0 via EIS data would need to have an overpotential value up to 50 mV. We thus suggest, before implementing EIS measurements at a fixed DC potential corresponding to the formal potential of the redox probe in a given electrochemical (electrode, electrolyte) system, the overpotential value to be evaluated by using Eq. 20 by running CV measurements at a medium scan rate (for example, 50 mV s −1 ).
Unlike A , which for a single electrode should be identical regardless of the method or the measuring conditions used, comparing each method using k 0 values of different redox probes in different electrolytes is not recommended since k 0 is a system dependent parameter, and is consequently inherently different for each redox probe/electrolyte system. Commenting on both approaches, CV-based methods provide more reliable results in comparison with the EIS-based method, where the validity of the results is dependent on the overpotential value of the examined system. Nonetheless, CV-based methods, while considered to be more reliable, they do need special attention especially in systems traversing between the quasi-reversible to the irreversible regime for all the reasons stated above.

Conclusions
We have applied and discussed the most widely used techniques and methods reported in literature for the calculation of electrode electroactive area ( A ) and heterogeneous electron transfer rate constants ( k 0 ) using graphite screen-printed electrodes. Considering the factors that should be taken into account when these methodologies are used for non-flat electrode surfaces [11,[13][14][15], as well as the approximation involved in graphical methods (for example, the Gileadi method), or the use of fitted values from impedance data, the reported results are considered as estimated values. We believe that we have prepared a useful and detailed guide which can instill the proper experimental demeanor in electrochemists.
The electrode electroactive area was calculated using both double-step chronocoulometry (Anson equation) and cyclic voltammetry (Randles-Ševčík equations). Once chronocoulometry with an outer-sphere redox probe (RuHex) produces symmetric and reproducible chronocoulograms and the respective Anson plots Q vs. t 1/2 are linear with an intercept equals (or it is almost) zero, its use is recommended. On the other hand, when cyclic voltammetry is employed, it is very important the use of the suitable form of the Randles-Ševčík equation with respect to the peak separation value(s) of the cyclic voltammograms at different scan rates. We also introduced " L ", a parameter dependent on the system's behavior (quasi-reversible or irreversible) for the correction of the experimentally obtained peak current values when A is calculated from the slope of a The heterogeneous electron transfer rate constants were calculated using both cyclic voltammetry and electrochemical impedance spectroscopy. Even though cyclic voltammetry-based methods (Nicholson, Klingler-Kochi, and Gileadi) are proved to be more elaborate, they can be applied to every system, considering the nΔE p for a given electrochemical system. As a rule of thumb, when nΔE p >150 mV the use of Klingler-Kochi method is suggested, while the Nicholson method remains the best choice when the nΔE p values are lower (ideally up to 105 mV). For an electrochemical system that its cyclic voltammetric behavior at increasing scan rates visibly transits from (quasi)reversible to irreversible regime, the graphical Gileadi method remains also a reliable alternative. On the other hand, electrochemical impedance spectroscopy can provide reliable values for electrochemical systems that exhibit low (typically < 50 mV) overpotentials and their impedance spectra can be modelled to a Randles circuit.
Funding Open access funding provided by HEAL-Link Greece.

Fig. 9
Current-overpotential curves generated by the (red line) Butler-Volmer equation and (black line) its linear approximation at an overpotential range of ± 10 to ± 300 mV. Highlighted areas designate the average current error percentage I error (%) ; green: 5%, yellow: 31%, orange: 69%, and salmon: > 86%. was set to 0.5 in the Butler-Volmer equation and i 0 was set to 10 −6 A in both cases